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High-throughput chromatin motion tracking in living 
yeast reveals the flexibility of the fiber throughout 
the genome 
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Chromosome dynamics are recognized to be intimately linked to genomic transactions, yet the physical principles 
governing spatial fluctuations of chromatin are still a matter of debate. Using high-throughput single-particle tracking, we 
recorded the movements of nine fluorescently labeled chromosome loci located on chromosomes III, IV, XII, and XIV of 
Saccharomyces cerevisiae over an extended temporal range spanning more than four orders of magnitude [10" 2 -10 3 sec). 
Spatial fluctuations appear to be characterized by an anomalous diffusive behavior, which is homogeneous in the time 
domain, for all sites analyzed. We show that this response is consistent with the Rouse polymer model, and we confirm the 
relevance of the model with Brownian dynamics simulations and the analysis of the statistical properties of the trajectories. 
Moreover, the analysis of the amplitude of fluctuations by the Rouse model shows that yeast chromatin is highly flexible, 
its persistence length being qualitatively estimated to <30 nm. Finally, we show that the Rouse model is also relevant to 
analyze chromosome motion in mutant cells depleted of proteins that bind to or assemble chromatin, and suggest that it 
provides a consistent framework to study chromatin dynamics. We discuss the implications of our findings for yeast 
genome architecture and for target search mechanisms in the nucleus. 



[Supplemental material is available for this article.] 

Chromosome higher-order architecture has been increasingly 
studied over the last decade thanks to technological breakthroughs 
in imaging and in molecular biology ranging from painting single 
chromosomes in fixed cells (Cremer and Cremer 2001) to mapping 
physical interactions between genomic elements at the genome- 
wide level (Hi-C) (Dekker 2008). It is now established that chro- 
mosome large-scale organization plays important roles in all as- 
pects of genomic transactions, allowing contact between distant 
chromatin loci (Gondor and Ohlsson 2009) that trigger concom- 
itant genomic processes in cis and in trans. Although several 
models have been proposed to describe the folding of chromo- 
somes (Lieberman-Aiden et al. 2009; Bohn and Heermann 2010; 
Bancaud et al. 2012), the physical principles governing this orga- 
nization are still largely debated. 

Despite its small size, budding yeast has become a unique 
model that recapitulates some of the main features of higher eu- 
karyotes, including heterochromatin-like regions, and an organi- 
zation mediated by the nuclear envelope (NE), the nucleolus, and 
the spindle pole body (SPB). Structural data obtained by statistical 
positioning of a gene in a yeast cell population led to a surprisingly 
simple model to define yeast nuclear architecture (Zimmer and Fabre 
2011): Chromosome position can be predicted by a few parameters 
such as genomic arm length, telomeres (TEL), and centromeres 
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(CEN) tethered to the NE via nuclear-envelope-tethered proteins 
and to the SPB via microtubules, respectively (Bystricky et al. 2005; 
Therizols et al. 2010; Zimmer and Fabre 201 1). This description was 
recently complemented by the first Hi-C comprehensive maps 
(Rodley et al. 2009; Duan et al. 2010), which confirmed an or- 
ganization guided by nuclear landmarks, including TEL that 
congregate in foci (Gotta et al. 1996; Schober et al. 2008). Recent 
Brownian dynamics (BD) simulations confirmed this structural 
model by recapitulating Hi-C and imaging data, assuming that 
physical tethering at TEL and CEN and volume exclusion were 
driving chromosome conformations (Tjong et al. 2012; Wong et al. 
2012). 

The large spatial fluctuations of chromosome loci inferred 
from particle tracking in living yeast suggested a highly dynamic 
character of chromosome architecture at the molecular level 
(Marshall et al. 1997; Heun et al. 2001; Hihara et al. 2012). Chro- 
matin dynamics appeared to be determined by nuclear constraints, 
in particular, the NE (Heun et al. 2001; Bystricky et al. 2004), and 
by the position along the chromosome, e.g., TEL (Bystricky et al. 
2005). Furthermore, a large body of data provides information on 
spatio-temporal dynamics in wild-type versus mutant cells (see, 
e.g., Bystricky et al. 2009) and elucidates how chromatin properties 
are regulated by structural proteins. However, the quantitative 
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analysis of these movements remains ChrlV 
controversial, and essentially two models ■ 
have been used: (1) Chromatin segments J 
were suggested to undergo normal Brow- 
nian fluctuations at small time scales 
and to be confined in regions of —0.3 
ixm 2 at time scales longer than —100 sec 
(Marshall et al. 1997); or (2) the move- 
ment of GAL 1-10 genes on chromosome 
II appeared to follow an anomalous be- 
havior characterized by a diffusion co- 
efficient of -0.4 (Cabal et al. 2006). This 
anomalous behavior was also detected in 
bacteria (Weber et al. 2010) and thor- 
oughly studied using polymer models. In 
yeast, the narrow temporal range of mea- 
surements within less than two decades was 
insufficient for a truly quantitative analysis 
of chromatin dynamics. 

Here we develop high-speed live 
cell imaging and high-throughput image 
analysis techniques to monitor the spatio- 
temporal fluctuations over an extended 
temporal range spanning more than four 
decades (10~ 2 -10 3 sec) of nine chromo- 
some loci on chromosomes III, IV, XII, and XIV of the yeast Sac- 
charomyces cerevisiae. We show that chromatin segments move 
subdiffusively and that the characteristics of the anomalous re- 
sponse are largely conserved for every locus distant by more than 
—50 kb from TEL and CEN. Similar dynamic behavior was also 
observed for two telomeric loci whenever they adopt a central lo- 
calization in the nucleus, thus leading to a consistent picture for 
chromosome movement within the entire nuclear volume. We 
then demonstrate that these dynamics are consistent with the 
Rouse model using BD simulations and statistical analysis of the 
trajectories and suggest that chromatin is a flexible polymer in 
yeast, characterized by a persistence length of <30 nm. Finally, we 
show that the Rouse model is also relevant in mutants for proteins 
involved in chromatin structure, and propose that this model de- 
fines a consistent framework to study chromosome motion. Our 
findings significantly advance the quantitative understanding of 
chromosome dynamics and their implications for yeast genome 
architecture. 
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Figure 1 . Collection of loci used to monitor chromosome spatio-temporal dynamics. Arrows indicate 
the genomic positions at which TetO or LacO repeats were inserted in four chromosomes. These loci are 
also represented by orange spheres on the yeast nuclear map obtained by Hi-C (Duan et al. 201 0). The 
two purple circles represent the telomeres that were tracked with the same labeling technique, and the 
green circles correspond to the centromere. The orange box corresponds to the rDNA locus, which 
represents —1 Mb of ribosomal DNA. The preferred localization of these loci was scored using a sub- 
division of the nucleus into two regions of equal surface (orange and green for central and peripheral 
localizations, respectively). 



2009). The motion of chromosome loci was recorded using a bright- 
field microscope and analyzed using a customized high-throughput 
tracking software based on the Multiple-Target Tracing algorithm 
(Supplemental Fig. SI; Serge et al. 2008; Albert et al. 2013). Because 
physical tethering of chromosomes may occur at the nuclear pe- 
riphery (Heun et al. 2001; Hediger et al. 2002), the nucleus was di- 
vided into two regions of equal surfaces, and every tracked locus was 
automatically assigned to a central or peripheral localization based 
on the segmentation of the first image of the acquisitions (see image 
in Fig. 1). 

We first focused on loci located in the nuclear center, far from 
constraints imposed by the nuclear membrane, and representing 
—60% of the observed positions (Fig. 2A). We then characterized 
their motion by plotting the temporal evolution of the mean 
square displacement (MSD), which describes the average dis- 
placement after a given time lag r: 



MSD(t) = <(r(t + r)- 



-r(t)) 2 >, 



(1) 



Results and Discussion 

The dynamics of chromosome loci are uniform throughout 
the nucleus and anomalous over a broad temporal range 

We assayed the dynamics of chromatin in a collection of yeast 
strains each bearing a single fluorescently labeled chromatin locus 
(see Methods for details). The labels were inserted on chromosomes 
III, IV, XII, and XIV at loci distant from TEL and CEN (Figs. 1, 2A). 
The chromosomes tested ranged in size from 300 kb to 1500 kb, 
which is well suited to assess whether the dynamics are affected by 
chromosome length. Motion was analyzed over a temporal range 
spanning more than four decades (15 msec to 400 sec) (Fig. 2B). This 
broad temporal analysis is more than two orders of magnitude larger 
than in previous studies performed in living yeast (Heun et al. 2001; 
Bystricky et al. 2004; Cabal et al. 2006) or in bacteria (Weber et al. 
2010), and investigations on similar time domains were only con- 
ducted in mammalian cells for telomere dynamics (Bronstein et al. 



with r the position vector, and (-> denoting the temporal average. 
We observed similar responses for every locus (data points in Fig. 
2B), which appeared to fold onto a single master curve (solid line). 
This curve follows an anomalous diffusive behavior over a broad 
temporal range characterized by a power-law scaling response in t a 
with a = 0.52 ± 0.08 (a = 1 for normal diffusion). This result is 
consistent with an earlier study describing the motion of GAL1-10 
genes, which exhibited a nonlinear MSD response characterized by 
an exponent of —0.4 in the temporal range of 4-80 sec (Cabal et al. 
2006), but it departs from other contributions showing normal 
diffusion and a plateau in the long time limit (see, e.g., Heun et al. 
2001). The anomalous subdiffusion response is associated with 
a slow increase in MSD over time. Its persistence over an extended 
period indicates that chromatin loci explore a broad region of 
the nucleus of ^/0.2 \±m 2 —450 nm in —400 sec. This estimate is 
consistent with the extent of "gene territories" of —500 nm, which 
are obtained by statistical mapping of the position of loci in cell 
populations of —2000 individuals (Berger et al. 2008). 
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Next we studied the motion of chro- 
mosome loci located at the periphery and 
observed that they frequently oscillated 
between the nuclear center and periphery 
during the time course of the acquisition 
(Supplemental Fig. SI), suggesting that the 
localization classification based on a single 
image is somewhat artificial. We did not 
detect significant differences in the MSD 
response for peripheral or central localiza- 
tion for the three loci on chromosome XII 
(see, e.g., position 240 kp in Fig. 2C) and 
on chromosomes IV and XIV (data not 
shown), so the anomalous diffusive re- 
sponse appears to be largely conserved and 
appropriate to describe chromosome seg- 
mental dynamics for loci distant by more 
than -50 kb from TEL and CEN. 
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Figure 2. Chromatin dynamics in living budding yeast. (A, upper panel) Fluorescence micrographs of 
different yeast strains for interframe intervals of 200 msec (for interframe intervals of 1 6-50 msec and 
mutants, see Supplemental Fig. S2); (lower plot) the spatial distribution of these loci in the nuclear 
volume. (B) The temporal evolution of the MSD for seven different chromatin loci with central locali- 
zation is plotted over a broad temporal range covering more than four orders of magnitude. These data 
show a universal behavior characterized by an anomalous diffusive response associated to power low 
scaling comprised between t 0AA and t 0 6 (thin dashed lines). (Black solid line) The Rouse regime, which is 
associated with an exponent of 0.5 (see details in Fig. 4). Note that each data set is the average over 30- 
200 trajectories and that we selected tracks with signal-to-noise ratios (SNR) >30 dB. (C) The MSD 
response for loci with central or peripheral localization is compared (yellow and green data sets, re- 
spectively), showing similar dynamics. (Black solid line) The behavior measured in B. The departure of 
the two curves in the short time regime is associated with the lower SNR of peripheral loci (data not 
shown). 



The dynamics of telomeres are 
anomalous in a central localization 

In yeast, the 32 TEL tend to cluster near 
the NE, and the motion of tagged chro- 
mosome sites in the vicinity of a TEL is 
constrained (Heun et al. 2001; Sage et al. 
2005). Using our reference behavior for 
chromosome loci located halfway along 
chromosome arms, we set out to perform 
an MSD analysis on TEL 3R and 14L (Fig. 1; 
hereafter denoted Tel3R and Tell4L). Al- 
though TEL predominantly accumulate 
near the periphery (Fig. 3 A, right panel), 
their positioning in the nuclear center 
can be detected with an occurrence of 
—30%, and we analyzed their dynamics 
in this central region. The amplitude of 
spatial fluctuations in the short time re- 
gime appeared to be relatively similar com- 
pared with the reference response (Fig. 3A), 
although Tel3R appeared to be slightly less 
mobile after —10 sec than Tell4L, pre- 
sumably because telomeric interactions 
between the two ends of chromosome 
XIV are much less frequent than that of 
chromosome III (Bystricky et al. 2005; 
Duan et al. 2010). 

In contrast to the bulk of chromo- 
somes, telomere dynamics appeared to be 
reduced at the nuclear periphery (Fig. 3B), 
in agreement with the idea that telomere 
peripheral localization is mediated by 
protein interactions involving, among 
others, Sir4 and Yku70/Yku80 (Hediger 
et al. 2002; Taddei et al. 2004; Bupp et al. 
2007; Schober et al. 2009). It has been 
shown that deletion of either YKU70 or 
SIR4 leads to derealization of some TEL 
to the nuclear center (Laroche et al. 1998; 
Hediger et al. 2002), but it remains un- 
clear whether this deletion affects TEL 
dynamics aside from chromosome III 
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chromosome segments dominate their 
local dynamics. This approximation, which 
leads to the "Rouse" regime (Rouse 1953), 
has been extensively documented for con- 
centrated solutions of synthetic polymers 
(see, e.g., Ewen and Richter 1997) and was 
recently validated in vivo for the bacterial 
chromosome (Weber et al. 2010). 

Chromosomes are modeled by a 
Rouse chain, which is composed of a se- 
ries of N elastic segments of length b. The 
motion of each bond is determined by its 
elasticity defined by the stiffness 3k B T/b z 
with k B T the thermal energy, and its vis- 
cous drag £. Assuming that chromosomes 
are tethered at CEN and TEL, the MSD of 
monomers can be solved analytically in 
the short time regime (see Model section): 
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12b z k B T 
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Figure 3. Telomere dynamics in living budding yeasts. (A, left plot) The MSD temporal evolution for 
telomere XIV-L and lll-R with central localization (red and blue data sets). (Black solid line) The average 
response obtained in Figure 2B. (Right p\ot) Telomeres preferentially accumulate at the periphery (green 
and orange bars, respectively). (B, left plot) The comparison of the MSD response for telomere XIV-L at 
a central or peripheral localization in wild-type cells (open yellow and green symbols, respectively) or in 
sir4 mutant strains (crossed symbols) favoring a central localization of telomere XIV-L. 



(Bystricky et al. 2009). We thus investigated the movements of 
Tell4L in wild-type and in SIR4 mutant strains. Deletion of SIR4 
leads to a predominant localization of Tell4L in the center (62%), 
likely due to loss of anchoring to the nuclear membrane. Moreover, 
telomere mobility is increased in SIR4 mutant strains when the 
tagged locus resides at the periphery and even more drastically at 
the nuclear center (Fig. 3B), confirming that physical interactions 
of the extremities of chromosomes are an essential player of their 
spatio-temporal dynamics. Overall, this study tends to support the 
proposition of uniform fluctuations throughout the nucleus for 
chromatin loci provided that they are sufficiently distant from 
anchoring regions. 

Chromosome segmental dynamics are consistent 
with the Rouse model 

To interpret our results with a polymer model, we make the assump- 
tions that chromosome conformation and dynamics are mainly 
determined by (1) the permanent tethering of CEN to the SPB and 
the more transient attachment of TEL to NE (Zimmer and Fabre 
2011), and (2) the level of crowding mediated by chromosomes 
and nuclear proteins. Crowding is expected to screen out hydro- 
dynamic effects, which otherwise induce long-range interactions 
in between chromosome segments (Doi and Edwards 1988), so 
that nearest-neighbor elastic interactions between consecutive 



The Rouse segmental dynamics are 
characterized by an anomalous diffusive 
response uniform along the chain, in good 
agreement with our data. We wished to 
test the predictions of the Rouse model 
further by extracting the histogram of the 
displacement after a time lag of 0.18, 1.0, 
and 1.6 sec from the trajectories of locus 
240 on chromosome XII (Fig. 4A). The 
histograms were first fitted with a normal 
diffusion model (Fig. 4A, left panel and 
equation in inset), which relies on one 
fitting parameter, the diffusion coefficient 
D. The diffusion coefficient of 4.8, 2.2, 
1.8 X 10" 3 |xm 2 /secwas very different for r = 0.1 8, 1.0, and 1.6 sec, 
respectively, showing the inconsistency of this model. We then 
used the formalism described in Guerin et al. (2012) to compute 
the histogram of the displacement for a Rouse chain tethered at its 
extremities (see expression in Model section). One single param- 
eter was sufficient to reproduce the three distributions (Fig. 4A, 
right panel), and the quality of the fitting was significantly im- 
proved (note the departure of the red fit from the conventional 
model), hence strongly supporting the relevance of the Rouse 
model to describe the segmental dynamics of yeast chromosomes. 

We then performed BD simulations of yeast chromosome 
dynamics in a volume the geometry of which is consistent with 
yeast nuclear size, and using the structure derived from Hi-C ex- 
periments as the initial state (Duan et al. 2010). We assumed that 
chromosomes were attached to the SPB at their CEN (i.e., neglecting 
the contribution of microtubules connecting SPB to kinetochores), 
and freely diffusing at their TEL (chromosome III is represented in 
the left panel of Fig. 4B). The computed MSD for loci distant from 
CEN (>400 nm; blue data set in the right panel of Fig. 4B) repro- 
duced the anomalous diffusive response associated with a power- 
law scaling of —0.54 over a temporal domain of four orders of 
magnitude that we observed experimentally (dashed line in Fig. 4B, 
right panel). The motion is more restricted for segments located at 
10-100 nm from the CEN (pink data set) versus at a distance com- 
prised between 200 and 400 nm (red data set). In addition, su- 
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Figure 4. Rouse model for chromosome segmental dynamics. (A) The two graphs represent the step 
distribution functions measured for three different time intervals, which are represented by dashed lines 
in Figure 2A. The genomic position of the locus on chromosome XII is 240 kb, and the inter-frame interval 
of the acquisition is 59 msec. Note the deviation of the distributions toward longer displacements for 
increasing time lags, which is consistent with the progressive exploration of the nuclear volume. The 
three curves are fitted to a normal diffusion model (equation in the inset) or an end-tethered Rouse 
model (Eqs. 4, 5) (left and right panels, respectively). (B) Brownian dynamics simulation was imple- 
mented with polymer chains and in a realistic yeast nucleus using the nuclear conformation obtained by 
Hi-C as initial state (Duan et al. 201 0). Chromosomes are tethered at their centromeres (green spheres) 
and free at their telomeres. The figure in the left panel represents the orthographic projection of 100 
conformations of chromosome III in the course of one simulation of 10 6 time points. The polymer is 
elongated, and its diameter is —300 ± 50 nm. The MSD is computed for loci localized near the SPB, or at 
200-400 nm and more than 400 nm from this anchor (purple, red, and blue data sets, respectively), 
showing that the response follows Rouse dynamics (black dashed line) provided that the locus is suf- 
ficiently distant from CEN. 



perposing the conformations of chromosome III during 10 6 time 
steps shows an elongated state of diameter —300 nm. This con- 
formation is associated with volume exclusion between chromo- 
somes (see also Supplemental Fig. S4B), an effect that was recently 
shown to be a key parameter of yeast nuclear architecture (Tjong 
et al. 2012). We thus conclude that the dynamics of chromosome 
loci in living yeast are accurately described by the Rouse model. 

The Rouse dynamics imply that chromosome segments 
search for target sites by compact exploration (de Gennes 1982), 
meaning that they systematically visit neighboring sites as they 
look for targets during repair or transcription activation (Condamin 
et al. 2007). Because the probability that chromosome loci find 
a nearby target is increased (Guigas and Weiss 2008), the Rouse 
dynamics are likely efficient to create specific chromosome contacts 
in cis or in trans over distances of less than —400 nm, which are 
explored in a few minutes (Fig. 2B). The notion of territoriality re- 
mains controversial, whether it is absent based on the observation 



that broken ends of chromosomes are free 
to search the entire genome for appropriate 
partners (Haber and Leung 1996), or partial 
as telomere colocalization and coordinated 
movement would indicate (Bystricky et al. 
2005). Future BD simulation studies should 
therefore evaluate whether and how com- 
pact exploration is related to chromosome 
territoriality in the yeast genome. 

MSD fluctuations suggest that 
chromatin is highly flexible 

The level of flexibility of a chromatin 
fiber is a matter of debate. In vitro single- 
molecule experiments have shown that 
the persistence length l p of nucleosome 
arrays is 20-30 nm (Cui and Bustamante 
2000; Bancaud et al. 2006; Celedon et al. 
2009; Kruithof et al. 2009). Conversely, 
analyses based on polymer models of in 
vivo distance measurements or Hi-C ex- 
periments (Bystricky et al. 2004; Dekker 
2008; Tjong et al. 2012) provided esti- 
mates of —200 nm and 66-134 nm, 
respectively. Physical modeling of the 
mechanical properties of nucleosome ar- 
rays predicts a persistence length smaller 
than 30 nm for chromatin (Ben-Haim et al. 
2001; Schiessel et al. 2001), except in the 
case of highly ordered arrays with strong 
nucleosome-nucleosome interactions 
(Wedemann and Langowski 2002), in 
which case l p increases to —250 nm. 

The amplitude of spatial fluctua- 
tions can be qualitatively linked to the 
flexibility of chromatin using the Rouse 
model, given that the bond length b and 
the viscous drag coefficient £ are equal to 
Zip and 6tvt\I p with r\ the viscosity, re- 
spectively (de Gennes 1979). Injecting 
this term in Equation 2 implies that the 
amplitude of the MSD increases with the 
square root of the polymer persistence 
length. The extent of spatio-temporal 
fluctuations therefore increases as the rigidity of the fiber increases, 
in qualitative opposition to free diffusion of an isolated tracer that 
increases as the diameter of the particle decreases. The MSD re- 
sponse can then be predicted for a persistence length comprised 
between 50 and 200 nm, given that the nuclear viscosity is r\ — 7 ± 
3 x 10~ 3 Pa.s based on GFP diffusion coefficient measurements 
(Kawai-Noma et al. 2006; Slaughter et al. 2007). This conjecture 
exceeds our measurements by a factor of 10-20 (dashed lines in Fig. 
5 A). In contrast, fitting of our data with the Rouse model leads to 
an amplitude of the MSD of 1.0 ± 0.2 X 10" 2 |jLm 2 /sec" 1/2 (solid 
line in Fig. 5 A), suggesting that the persistence length of chromatin 
in living yeast is —5 nm. This length scale is coherent with the 
dimensions of a single nucleosome, which is a — 5-nm by — 10-nm 
disk-like structure. Hence, the analysis of chromosome motion 
suggests that chromatin is a flexible polymer in vivo. 

Because the Rouse model is a simplified model that does 
not take volume exclusion between monomers into account, 
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Figure 5. Assaying chromatin flexibility from MSD analysis. (A, gray) The MSD response plotted in 
Figure 2B is analyzed with the Rouse model for different values of the persistence length of 200, 50, and 
5 nm, as specified in the inset. The dynamics of chromosome loci are consistent with highly flexible 
chromatin. (B) Brownian dynamics simulations of an isolated chain containing 100 segments have been 
performed to confirm the relevance of the Rouse model. Dimensionless variables are defined using l p and 



jectories of a locus on chromosome XII 
(position 680, Fig. 5C). Because the 
difference between the average and 
the maximal or minimal MSD is about 
threefold, this variability is unlikely re- 
lated to changes in persistence length, 
which should vary over two orders of 
magnitude from 0.3 to 30 nm. We rather 
suggest that this variability is associated 
with the existence of dynamic chromo- 
somal contacts, as previously indicated in, 
e.g., the dynamic random loop model de- 
rived to interpret distance measurements 
in human cells (Mateos-Langerak et al. 
2009). Dynamic and random contacts 
should transiently restrict spatial explo- 
ration (Supplemental Fig. S5), reduce the 
effective diffusion of a genomic site over a 
population average, and hence lead to an 
apparent persistence length smaller than 
chromatin intrinsic properties. We pro- 
pose that chromosome loci with maximal 
mean square displacements in the short 
time regime represent the "interaction- 
free" limit, because their motion is not 
hindered by distant contacts. Given that 
the difference between the average and 
the maximal MSD is less than threefold 
(Fig. 5C), we can evaluate an upper limit 
for the persistence length of 3 2 X 3 nm 
<30 nm. Consequently, chromatin ap- 
pears as a flexible polymer, and its per- 
sistence length can be qualitatively esti- 
mated to be less than —30 nm. 



r]l p h- k B T as characteristic length and time scales; the MSD of a central monomer is plotted with dif- 
ferent persistence lengths of 1 5, 30, 60, or 90 nm and for two different viscosities of 1 0" 3 and 1 0 2 P 
(blue data sets). These data sets fall on a master curve, which follows the power law scaling of the 
analytical Rouse model: The black dashed line is derived from Equation 2 using the dimensionless 
variables (see exact expression in inset). (C) The histogram of the average square displacement after 
0.1 , 0.6, and 4.0 sec is represented for —70 trajectories of a locus at position 680 on chromosome XII, 
with a broad cell-to-cell variability within a single clone. The lines are Gaussian fits and serve as guides 
to the eye. 



we wished to strengthen our conclusions by running BD of 
a single chain composed of self-avoiding segments (see Model 
section). The extracted MSD shows an anomalous diffusive 
response, as well as the dependence on the persistence length 
and the viscous drag (Eq. 2) expected from the Rouse model (Fig. 
5B). This analysis showed that the numerical pref actor of Equation 
2 is underestimated by 1.4, suggesting that volume exclusion 
enhances the rate of spatial exploration. This difference implies 
that the persistence length of —5 nm is overestimated, and 
a correct estimate may be as small as ~3 nm. This length scale, 
which is comparable to the DNA diameter, is unrealistic and 
points to an additional parameter related to the chromosomal 
environment. 

Each MSD curve represented in Figure 2B is the average re- 
sponse over a population of —100 cells. The clonal variability of 
the spatio-temporal dynamics is not shown in this graph, although 
this variability is large, as inferred from the distribution of average 
displacements after a time lag of 0.1, 0.6, and 4.0 sec for —70 tra- 



Chromatin structural proteins alter 
chromosome mobility 

Finally, we investigated whether the dy- 
namics of chromosomes was altered in 
the absence of chromatin structural pro- 
teins or chromatin assembly proteins. We 
first considered the HMG-like protein 
Spt2 in S. cerevisiae, which interacts with 
the histones H3 and H4 (Horn et al. 2002). The statics and the dy- 
namics of one locus on chromosome XIV were assayed in wild-type 
(WT) and mutant spt2 cells, showing an accumulation at the pe- 
riphery for the mutant and similar MSD responses in the short time 
regime (Fig. 6, right and left panels, respectively). We then per- 
formed the same analysis with Hholp (Levy et al. 2008), which is 
considered as the linker histone analog in S. cerevisiae, although its 
structure differs from the linker histone in other organisms be- 
cause it bears two globular domains. So far it remains unclear 
whether and how this protein participates in chromatin conden- 
sation. The deletion of HHOl increases the proportion of spots 
detected in the nuclear center compared with wild type, and in- 
creases — 1.5 -fold the amplitude of chromatin loci displacements 
(Fig. 6). Given that the ratio Hholp versus nucleosome remains 
poorly estimated (comprised between 1:37 and 1:4) (Freidkin and 
Katcoff 2001; Downs et al. 2003), these enhanced dynamics may 
be associated with the reduced frequency of chromosome in- 
teractions mediated by the two domains of Hholp and/ or with an 
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Figure 6. Chromatin structural/modifying proteins and chromosome dynamics. (Left plot) The MSD 
temporal evolution for a locus on chromosome XIV, genomic position 240 kb. (Black solid line) The wild- 
type situation was identical to the Rouse response described in Figure 3A. The dynamic response 
obtained in hhol (red symbols), asfl (blue symbols), and spt2 (purple symbols) mutant yeast. The two 
dashed lines represent fits to the data with the Rouse model. (Right plot) Asfl p and Hhol p participate in 
the nuclear localization of the tagged locus (the yellow and green bars correspond to central and pe- 
ripheral localizations, respectively). 



The folding of interphase chromo- 
somes in metazoans involves loops of var- 
iable sizes defining topological domains of 
0.2-5 Mb (Nora et al. 2012; Sexton et al. 
2012). Because the dimension of these 
domains, which can be described as 
end-tethered chromosome fragments, 
is consistent with the length of yeast 
chromosome arms, we speculate that the 
Rouse model is relevant to higher-order 
eukaryotes. This model has recently been 
implemented in numerical simulations 
(Bohn and Heermann 2010), showing that 
the dynamics of loop formation may in- 
terfere with the segmental motion of hu- 
man chromosomes, a hypothesis that re- 
mains to be tested quantitatively in vivo. 

Model 



onset in flexibility due to conformational changes at the nucle- 
osome level. Finally, we probed the role of the nucleosome 
assembly chaperone Asflp, which is involved in replication- 
dependent and replication-independent chromatin assembly 
(Loyola and Almouzni 2004). As for hhol deletion mutant, the 
labeled locus appeared to be more frequently central than in wild- 
type nuclei, and its mobility was 2.2-fold increased. The Rouse 
model therefore seems to account for the dynamics of chromo- 
some loci in WT cells, as well as in hhol, asfl, and spt2 mutants, 
suggesting that this approach is sufficiently generic to perform 
systematic motion analyses. 

Conclusion 

We demonstrate that the motion of chromosome loci in living 
yeast is characterized by anomalous subdiffusion with an expo- 
nent of 0.52 over a temporal domain spanning more than four 
orders of magnitude. We then show that this response, which 
appears to be relevant for the analysis of WT and mutant cells, is 
consistent with the Rouse dynamics. Because the Rouse regime is 
characteristic of polymers in crowded environments, our result 
strengthens the relevance of polymer models for the description 
of the yeast genome (Gehlen et al. 2012; Tjong et al. 2012; Wong 
et al. 2012). In addition, the Rouse dynamics is associated 
with compact exploration, implying that chromosome loci ef- 
ficiently search for nearby targets. We speculate that this prop- 
erty guides chromosome interactions in cis and trans, and 
envision that this hypothesis can be tested using BD simulations 
in combination with Hi-C experiments aiming to evaluate the 
number of chromosome contacts as a function of the genomic 
distance. Furthermore, the amplitude of chromatin spatio-temporal 
fluctuations shows that this polymer is flexible. Yet the cell-to- 
cell variability of the amplitude of chromosome movements 
only allows for a qualitative determination of the persistence 
length of <30 nm. We propose that this variability is associated 
with dynamic interactions between chromosomes, and future 
BD simulations may confirm or invalidate this conjecture, in 
turn, providing new insights on the statics and the dynamics of 
the yeast genome. 



Dynamics 

The motion of an end-tethered Rouse chain 
can be studied analytically using normal 
mode decomposition (Vandoolaeghe and Terentjev 2007), and the 
MSD of the nth segment (MSD n ) reads: 



MSD„(t) 
4Nb 2 



"7r 



P \sin 



Tt)' 



(3) 



with 



p the relaxation time of each mode j p - 



T R /p 2 and = L\l 



3tt k B T the Rouse time identified as the longest relaxation time of 
the chain. Note that this equation is valid in three dimensions 
(3D), and its two-dimensional (2D) version is readily obtained with 
a multiplication by 2/3, assuming that the motion is isotropic. The 
temporal evolution of the MSD is plotted in Supplemental Figure 
S4As. In the short time limit (t < tr), Equation 3 is dominated by 
the terms with large p, and it can be simplified to obtain Equation 
2. In addition, the formalism described in Guerin et al. (2012) (see 
details in the Supplemental Material) enables us to extract the 
distribution of the displacement R for the nth segment after a time 
lag t: 



Pn(R,r) = 



*n(r)=^SLsin* [p 



*n(r) 



exp- 



R z 
2^ n (r) 



(4) 



N 



N+lJ 



1- exp 



2£b 2 r 
' k B T 



1- cos [p 



N+lJ 



1 - cos p 



N+l 



(5) 



The movements of chromosomes, which were modeled by 
a series of bead linked by springs or ball-in-socket joints, was also 
analyzed by BD. The size of each bead was set to 15 nm in di- 
ameter, and it was assumed to contain approximately three 
nucleosomes, so that the full genome consists of 27,000 beads. 
Chromosomes were constrained in a sphere of 2 |xm in diameter, 
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and the CEN were fixed to the SPB using harmonic constraints. 
We used the Langevin dynamics algorithm implemented in the 
NAMD software (Phillips et al. 2005) with an appropriate in- 
house force field. The Langevin equation for the nth monomer 
has the form: 

m >< -fit (f) = ~Srad x E + f„(t)-m n b n -£ (t ), (6) 

where the random force f n (t) is derived from a Gaussian distribu- 
tion with the properties: 

(fn(t)) = 0 (7) 



Cell culture 

Cells were grown in YPD or YNB at 30°C starting from a fresh 
overnight culture. They were diluted at 10 6 cells/mL, and harvested 
when OD 60 o reached 4 X 10 6 cells/mL and rinsed twice with the 
corresponding SC media. Cells were then spread on a SC media 
patch containing a 2% agarose and 2% glucose coated slide. A 
coverslip was deposited over the cells, and the construction was 
finally sealed with "VaLaP" (one-third vaseline, one-third lanoline, 
one-third paraffin). Live microscopy was limited to 20 min after 
mounting the coverslip. Extreme care was taken to ensure that cells 
are at approximately equal stages of growth based on their shapes, as 
inferred from a single image of transmission light microscopy, and 
based on the nucleus shape using the roundness parameter i?, which 
was defined as: 



(f n (t)f n (0)) = 2k B T{m r fi(t), 



(8) 



with £ the friction coefficient. The energy £ of each bead is the sum 
of three contributions associated with the confinement in the 
nucleus, the elasticity of the chain, and the repulsive interactions 
between monomers, which are expressed as: 



■ + E V 



Eamf{r<ri) = 0 and E conf (r>ri) = E x (r-rif, 



(9) 



(10) 



with r the distance of the bead to the center of the nucleus, and r ± 
the radius of the nucleus. 



E s pring(n, n+l) = k s {r n -r n+l ) 



Evdwihj) = U 



(ii) 



(12) 



The parameters of this model can be found in Table 1 . BD were 
run during 3 X 10 7 time steps, starting from the conformation 
obtained by Hi-C (Duan et al. 2010) equilibrated during 10 6 time 
steps. We also performed BD to investigate whether spatio-temporal 
fluctuations were different for 1, 4, or 10 segments, or for one seg- 
ment of large diameter along the chain. We did not detect differ- 
ences on the MSD response (data not shown), in agreement with our 
observation of similar dynamics with Lac or Tet labeling, which 
measure —10 and 7 kb in length, respectively. 



Methods 

Plasmids and strains 

All the strains we used are described in Supplemental Table 1. 

Table 1 . Set of parameters to run Brownian dynamics simulations 
T 300 K 



u 

m n 
b n 



1 kcal mol nm 

7 nm 
-0.01 kcal mol 1 
10 a. u. 
10 a. u. 
0.01 kcal mol 1 nm 



R = 477 X Surface/ Perimeter 2 



(13) 



with Surface the number of pixels after segmentation of the 
nucleus, and Perimeter the outline in pixels (see Supplemental 
Fig. S2). We considered nuclei with R > 0.8 to avoid mitotic 
phenotypes. 

Live cell imaging 

Yeast cells were imaged at room temperature with a BX-51 upright 
microscope equipped with a laser diode (Lumencor), a 100X oil 
immersion objective (NA = 1.4), and an EMCCD camera (Andor 
DU-897), as described in Hajjoul et al. (2009) (see Supplemental 
Video 1). The excitation emission at 470 ± 10 nm was set to 7.53 
W/|xm 2 , about 100 times less than the value that is necessary for 
normal cell growth (Carlton et al. 2010), and we observed con- 
secutive rounds of divisions during ~2 h (data not shown). For 
slow acquisitions of interframe intervals of 2 sec or more, we used 
a 100X oil immersion objective with a variable NA set to ~ 1.1-1.2. 
The acquisition was performed with a shutter synchronized with 
the camera to minimize photobleaching, and we checked for drifts 
by monitoring the position of the center of the nucleus in the first 
and last image of the recording. Optimal interframe intervals of 
—20 msec were reached by cropping regions of interest and using 
a 2 * 2 binning with the Andor iQ imaging software. The loss in 
pixel size resolution was compensated with a 2x lens placed in the 
microscope light path. 

We focused on 2D single-particle tracking and checked that 
similar information was retrieved with 3D tracking using micro- 
mirrors (Supplemental Fig. S3; Hajjoul et al. 2009, 2011). We note 
that 2D acquisitions may bias MSD responses for loci preferen- 
tially localized at the periphery, because in some cases central loci 
in the projected view may be peripheral in the 3D nucleus. 

Data analysis 

Video sequences consisted of 300 consecutive images, and we dis- 
played MSD traces on 150 time intervals at most in order to ensure 
the statistical relevance of mean displacements. The movements of 
chromosome loci were systematically analyzed using a custom 
software that was developed in Matlab (Mathworks) (Serge et al. 
2008; Albert et al. 2013 [source code at http://jcb.rupress.org/ 
content/202/2/201/suppl/DCl, and the executable file at ftp:// 
intermtt:MTTinterface@ftp.laas.fr/]). This software enabled us to 
extract (x, y) coordinates by Gaussian fitting, to reconstruct the 
trajectories, and to compute the MSD and the step distribution 
function (see Supplemental Video 1). Note that the MSD in fixed 
cells was 20- to 30-fold lower than in living cells (Supplemental 
Fig. S3B). 
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